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Abstract. We present VLT-ISAAC L-band spectroscopy toward 10 stars in SVS 4, a 30" x 45" dense cluster of pre-main 
sequence stars deeply embedded in the Serpens star forming cloud. The ISAAC spectra are combined with archival imaging 
from UKIRT and ISOCAM to derive accurate extinctions toward the SVS 4 stars. The data are then used to construct a spatial 
map of the distribution of ice in front of the cluster stars with an average angular resolution of 6" or 1500 AU, three orders 
of magnitude better than previous maps. We show that water ice is present throughout the region and confirm the presence of 
methanol ice with an abundance of up to 25% relative to water. It is shown that methanol ice maintains a very high abundance 
relative to PL throughout SVS 4, but drops by at least an order of magnitude only 75" away from SVS 4. The maps indicate that 
some of the lines of sight toward the SVS 4 stars pass through the outer envelope of the class protostar SMM 4. The abundance 
of water ice relative to the refractory dust component shows a sudden increase by 90% to (1.7 ± 0.2) x 10~ 4 relative to H 2 at a 
distance of 5 000 AU to the center of SMM 4. The water ice abundance outside the jump remains constant at (9 ± 1) x 10 5 . 
We suggest that this is an indication of a significantly enhanced ice formation efficiency in the envelopes of protostars. The 
depletion of volatile molecules in the envelope of SMM 4 is discussed. In particular, it is found that up to 2/3 of the depleted 
CO is converted into CO2 and CH3OH in the ice. Therefore, only 1/3 of the CO originally frozen out will return to the gas phase 
as CO upon warmup. 

Key words. Astrochemistry, Stars: formation, circumstellar matter, infrared: ISM, ISM: molecules, dust, extinction 
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03 1 . Introduction infrared source, most of these observational studies have con- 

' centrated on ices along single, isolated lines of sight. Large 

.£ In the last decade, it has become increasingly clear that the variations in the abundances f some ice species have been 

^ . interaction between molecular gas and ice mantles on dust obseryed ^ diflferent regionSi ^ ^ lack of informa _ 

H ; grains in dense clouds significantly affects the overall chem- tion hag often made any interpretation difficult Examples of ice 

LU , istry of the cloud as well as the way molecules are observed. . , , ■ , . . ■ 

3 1 species with strongly varying abundances relative to water ice 

The chemistry of a dense cloud is influenced by the pres- , ■ , • \ T J . . miinnrl fn I ~a TTI 

1 include methanol ice (Dartois et al. 1999; PontorJDidan et al. 

ence of ices through reactions taking place on the surfaces of 2003a)> QCN" ice ( Pendleton et al. 1999; Gibb et al. 2004) and 

dust grains and by changing the composition of the gas when CQ ice ( Pontorjrjidan et al 2003b , Qther abundant ice species 

ice species are desorbed through heating or energetic process- m found tQ haye a f air iy constant ab undance relative to wa- 

ing < Aikawaetal. || 2003 | ) . Additionally, molecules may be re- fcr ^ such afj ^ (G erakines^5]S. This shows that 

moved from the gas-phase by freezing out onto dust grains i n ^ composition of interste llar ice depends on the physical con- 

the early phases of star formation (e.g. | j0rgensen et aH | 2002|) . ditions under which ft formed ^ fof a complete smdy of 

Many observational studies have attested to the high abun- ^ chemical evolution of interst ellar ice, spatial information 

dance of a large range of molecular ice speci es in dark clouds ag weU ag information on the physical environment in 

and in star forming regions in particular (e.g . Gibb et al. 2004; which ft h found ^ needed Infrared spectrometers are now 

Booaert & Ehrenfreund 2004; van Dishoeck et al. 2003). Since ... , , ■ . , ,. f ■ , . . ■ 

1 6 * " sensitive enough to record spectra along lines of sight toward 

interstellar ice can only be observed in absorption toward an,, a ■ r a u- u u * 

J F closely spaced infrared sources, which can be young stars or 

~Z ~, „. . , ~~~ „ t,. „ . ., .. background stars. This opens the possibility of directly obtain- 

Sencl offprint requests to: K. M. Pontoppidan e-mail: r r J J 

pontoppi@strw. leidenuniv.nl in S s P atial information on the distribution of interstellar ices. 

* Based on observations obtained at the European Southern For example, |Murakawa et alj (|2000j) constructed a water ice 

Observatory, Paranal, Chile, within the observing program 71.C- map of the Taurus molecular cloud using bright background 

0252(A). stars. The resulting spatial resolution of this study was ~1 line 
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of sight per 100 square arcminutes. In this article we present 
high angular resolution (1 line of sight per 0.02 square arcmin- 
utes) observations of the distribution of ices in a small region 
of the Serpens star forming cloud centered on the cluster SVS 
4. 

SVS 4 is a small but dense cluster of low- to intermediate 
mass (L - 1-50 L ) pre-main sequence stars located in the 
south-east ern core of the Serpens molecular cloud. It was first 
studied bv lEiroa & Casalil {l989), who found it to be one of 
the densest YSO clusters known with 1 1 stars within a region 
only 12000AU across, corresponding to a stellar density of 
~ 5x 10 4 M pc" 3 . One of the cluster members, SVS 4-9/EC 95, 
is the brightest X-ray source in Serpens and is probably a very 
young int ermediate mass star with abnormal coronal activity 
(Preibisch 2003). The center of the cluster is located only 30" 
or 750 AU from the deeply e mbedded class protostar, SMM 
4 (e.g. Hogerheiid e et alJl999t) . Lines of sight toward the stars 
in the SVS 4 cluster may therefore intersect part of the outer en- 
velope of SMM 4. This system thus provides a unique oppor- 
tunity to directly probe the amount of freeze-out near deeply 
embedded young stella r objects. Detailed millimeter studies 
of pre-stellar cores fe.g. lCaselli et alJl999HBergin et al| 2002: 
iTafalla et alj E004) and protostars (e.g. lj0rg^nsenetal] l2OO4) 
indicate that a significant fraction of the gas-phase molecules 
(> 90% of the condensible species) is depleted onto grains in 
the cold parts of the envelopes where the timescales for freeze- 
out are shorter than the age of the core. For the case of proto- 
stellar envelopes, J0rgensen et al. propose a 'drop' abundance 
structure where the gas-phase abundances are relatively high in 
the low density outermost envelope and then drop abruptly to 
a very low value at the radius where freeze-out becomes effec- 
tive. In the inner warm regions of the envelope, the ices evapo- 
rate resulting in high gas-phase abundances. Our observational 
ice maps can test such abundance structures directly. 

Pontoppi dan et al.l J2003al) presented low-resolution L- 
band spectra toward two of the SVS 4 cluster members, SVS 
4-5 and SVS 4-9, and found that the ices in front of these 
two stars contain a very high abundance of methanol ice of 
up to 25% compared to water. This was the first detection of 
methanol ice outside of high-mass star forming regions and is 
among the highest abundances observed. Previous searches for 
methanol ice in low-mass star forming regions have only pro- 
vided strict upper limits to the abun dance of meth anol down 
to a few % relative to water ice JChiar et alj fl9961 The rea- 
son for this large variation in the methanol ice abundance is 
currently not understood. Recent laboratory experiments have 
shown that methanol can be formed in dense clouds through 
successive hydrogenati on of CO a t 10-20K on the surfaces of 
dust grains dWatanabe et alJl2003l) . It is found that the forma- 
tion of methanol ice along this route depends primarily on the 
abundance of atomic hydrogen, the thickness of the ice mantle 
and the temperature of the ice, with the most efficient forma- 
tion of methanol at 15 K. Thus, a variation in these parameters 
could be related to the observed differences in methanol ice 
abundance. Alternatively, bombardment with energetic protons 
has also been found to produce moderate yields of methanol in 
CO:H 2 ice mixtures (Hudson & Moorel ll999). The study of 
the relation of interstellar methanol ice to the cloud environ- 



ment in which it is found therefore provides direct constraints 
on the formation mechanism of this important molecule. Due 
to the high stellar density, the proximity to SMM 4 and the very 
high methanol ice abundance of the SVS 4 cluster, it presents 
an excellent case for studying the distribution of methanol-rich 
interstellar ice on small angular scales in a protostellar environ- 
ments. 

In this paper we present low-resolution L-band spectra of 
10 SVS 4 stars obtained with the ISAAC spectrometer on the 
Very Large Telescope (VLT) at Paranal in Chile. Additionally, 
we use archival data from the Infrared Space Observatory (ISO) 
and the United Kingdom Infrared Telescope (UKIRT) to fur- 
ther constrain the environment. We use the spectra to produce 
a detailed map of the distribution of water and methanol ice in 
the SVS 4 region with a spatial resolution of 6"=1500 AU. The 
corresponding line of sight density is one per 0.02 square ar- 
cminutes, which is more than three orders of magnitude higher 
than previous ice maps. 

This paper is organised as follows: §|2]describes the obser- 
vations from ISAAC as well as archival data used. §0 shows 
how the basic observational parameters have been derived. §|4] 
interprets the relation to the SMM 4 envelope and §[5]discusses 
the distribution of ices in the observed region. 

2. Observations and archival data 

L-band spectra of 10 sources in the SVS 4 cluster were ob- 
tained using the ISAAC spectrometer mounted on UT 1 of the 
VLT on the nights of July 8-9, 2003. The low resolution grat- 
ing and the 076 slit were used, yielding a resolving power of 
A/AA ~ 600 and an instantaneous wavelength coverage from 
2.84^4. 15/vm. Typical on-source integration times were ~40 
minutes per pointing. Additionally, the two brightest sources, 
SVS 4-5 and SVS 4-9, were observed in one setting centered 
on the CH-stretching mode of solid methanol at 3.53 fim us- 
ing the medium resolution grating with a resolving power of 
A/AA ~ 3300. The spectra were reduced using IDL fo llow- 
ing the procedure described in Pontop pidan et alJ (2003bJ. The 
spectra were ratioed by the early-type standard stars, BS 6629 
(A0V) and BS 7236 (B9V) to remove telluric absorption lines. 

A 4.07 yum image was obtained by co-adding 8 acquisition 
images from the spectroscopic observations and was flux cali- 
brated relative to the acquisition image of the standard star BS 
7348 (B8V). The 4/mi flux of the standard star was estimated 
by extrapolating the V-band magnitude of 3.95 to 4 /urn using 
blackbody colours, yielding M\ lim = 4.15. The spectra were 
then scaled to the 4.07 //m photometric points. We estimate the 
accuracy of the photometry to be better than 0.1 mag. Standard 
stars of spectral type A or B are known to have photospheric 
hydrogen lines in absorption. This may result in excess line 
emission in the ratioed spectra. The two standard stars show 
the presence of the Bra line at the 7% level relative to the con- 
tinuum while other hydrogen lines have strengths of < 5%. The 
Bra line was not detected in any of the SVS 4 spectra and any 
residual emission lines are due to the standard stars. Since the 
hydrogen lines do not affect the ice bands, no attempt to correct 
for them was made. 
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J- H and ZT-band imaging from UFTI at the United 
Kingdom Infrared Telescope (UKIRT) were extracted from the 
UKJRT archive, and ISOCAM 6.7 and 14.3 //m imaging as well 
as a CVF field centered on SVS 4 were extracted from the ISO 
archive. The near-infrared images were flux calibrated using 
bright 2MASS stars in the field. Standard aperture photometry 
was performed using 1.5" apertures. The absolute photome- 
try was found to be better than 0.1 mag when comparing to 
2MASS. The ISOCAM data were reduced using the CIA re- 
duction package version DEC01. Because the SVS 4 sources 
are partially blended at the ISOCAM resolution, the raster 
maps were resampled to a finer sampling by a factor of two. 
ISOCAM fluxes were then derived from the broad-band images 
by fitting a PSF derived from an isolated source in the field to 
source positions taken from the 4.07 fim ISAAC image. The 
ISOCAM images where then flux calibrated using conversion 
factors of 2.32 adu g -1 s _1 mJy " 1 and 1.96adug~' s" 1 m Jy 1 for 
LW2 and LW3, respectively ( Blom maert et alJ boOOl) . where 
adu is the Analog to Digital Unit, and g is the gain. Any 
colour corrections to the broad band ISOCAM fluxes were 
found to be less than 2%. CVF spectra were extracted using 
3 pixels x 3 pixels = 18" x 18" apertures. 



3. Derivation of observed physical parameters 

3.1. Distance to the Serpens molecular cloud 

There is some controversy in the literature regarding the 
distance to the Serpens cloud. Distance estimates based 
on accurate phot ometry of a few stars (~10) range from 
250 nc JChavarria-K.et al1ll988h to 400 pc dChiatl Il997t 
iHog erheiide et al. 1999f). However, accurate photometric sur- 
veys of larger samples of stars seem to conv erge on the smaller 
distance of ~ 250 pc. Straiz vs et alJ {l996) find a distance of 
259 ± 37 pc using 105 stars. This result is supported by a pre- 
liminary distance estimate of 220-270pc using 2MASS pho- 
tometry of a similarly sized sample (Knude et al., in prep). In 
this paper all distance-dependent quantities are scaled to a dis- 
tance of 250 pc. 

3.2. Determination of extinction and water band 
continuum 

In order to study the distribution of ices toward the SVS 4 clus- 
ter, accurate optical depths of the water bands as well as col- 
umn densities of the refractory dust component must be de- 
termined. The optical depths of the water bands depend sensi- 
tively on the continuum chosen. However, since the blue wing 
of the water band is outside of the atmospheric L-band win- 
dow, the determination of a continuum will always have some 
degree of uncertainty. One solution is t o use a A"-band spectrum 
to define an empirical continuum (e.g. lMurakawa etlril EoOO). 
Unfortunately, only photometric JHK points are available for 
the SVS 4 stars. However, as will be shown in this section, 
most of the sources in SVS 4 have almost no infrared excess. 
This allows reasonably accurate modeling of the continua us- 
ing Kurucz stellar atmosphere models. Furthermore, a simulta- 
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Fig. 1. 7, H, 4 /mi colour composite o f SVS 4. The co ntours 
show the 850/mi SCUBA map bv bavis et ail Jl999l) . The 
cross indicates the central position of SMM 4. 

neous determination of the extinction toward each source can 
be derived with confidence. 

Since the effective temperatures of low-mass pr e-main se- 
quenc e stars are quite insensitive to their masses dSiess et alJ 
2000), Kurucz models with a low logg and T e g = 3500 - 
4500 K can be used as a rough model of the photospheric emis- 
sion. In the case of SVS 4-9 and SVS 4-10, the effe ctive tem- 
peratu res have been determined spectroscopically by Preibisch 
( 1999). Any infrared excess is modeled by a single blackbody 
of T = 500 - 1500 K or by a disk model from Dulle mond et alJ 
if the excess is large. The observed SED of each source 
was de-reddened to fit a Kurucz model at the J, H and K points 
and a blackbody or disk model was added to fit any infrared 
excess at longer wavelengths as well as the 3.8-4. 15/im slope 
of the L-band spectra. In cases where a strong infrared excess 
would also affect the near-infrared points, the extinction and 
infrared excess contribution were iteratively corrected until a 
fit was found, which reproduced both the photometric points 
and the slope of the L band spectra. This approach was chosen 
to provide a realistic continuum at 3 pm rather than to yield an 
unambiguous excess model. The fits to sources with a strong 
excess are therefore degenerate to some extent. However, due 
to the very high absolute values for the extinction (Aj > 5) 
the determination of water ice optical depths and extinction 
was found to be very insensitive within the allowed parame- 
ter space, and the uncertainties given on the derived quantities 
reflect the fitting degeneracy. 

The extinction law at near-infrared wavelengths is well 
described by a power law, Aj/Aj - (A/ 1.25 um) a , wher e 
a = 1.6 - 2.0 dCardelli et alJll989t iMartin & Whittetlll990h . 
A high value of a has been found t o best describe dust in 
some dense molecular clouds (e.g. Martin & Whittet 1990; 
iLutz et alJfl996l) . A low value for the power-law index of 1.6 
was attempted, but this was found to give extinction-corrected 
2-5 /mi fluxes, which were too steep to correspond to a black- 
body or Kurucz model for stars with no infrared excess when 
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Fig. 2. The SVS 4 /, //, 4/mi colour composite overlaid with 
ISOCAM 6.7 /an contours. 
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Fig. 3. The SVS 4 J,H,4/j.m colour composite overlaid with 
ISOCAM 14.7 fim contours. 

the {J - H) colours were used to determine the extinction. 
In particular, an extinction law index of at least 1.85 was re- 
quired to simultaneously fit all spectral and photometric points 
of sources like SVS 4-3, SVS 4-5 and SVS 4-9. A higher 
index could also be used since this would simply create a 
small infrared excess. Here, we adopt the minimum value of 
a - 1.85. In the mid-infrared ran ge, the emp i rical e xtinction 
law toward the galactic center by iLutz et al ] ill 9961) is used. 
This is consistent with the ISOCAM extinction law for Ser pens 
(At w? = 0A1A K and A LW 3 = 0.36A^) derived bv Kaas et alJ 
( 2004). The resulting SEDs and water bands are shown in Figs. 
@]|5]and|6] The corresponding fitting parameters are shown in 
Table |2 The extinctions are given in terms of the extinction 
in the /-band, A y, to avoid making unnecessary assumptions 
about shape of the extinction law in the optical wavebands by 
extrapolating to Ay. 



It is seen that the extinction toward SVS 4 is in general 
very high (Aj = 5 - 25), roughly corresponding an Ay of up 
to 100 mag. Most of the sources exhibit very little infrared ex- 
cess shortwards of 6 fim with the exceptions of SVS 4-6 and 
SVS 4-10. This results in a fairly accurate absolute determi- 
nation of Ay with an uncertainty dominated by the chosen ex- 
tinction law rather than the accuracy of the photometry or, in 
most cases, the SED model assumptions. We estimate a rela- 
tive uncertainty on Ay of 5%. A different extinction law may 
cause systematic effects larger than this. The bolometric lumi- 
nosities of the stars can then be estimated by integrating over 
the fitted Kurucz models. The derived stellar parameters are in- 
cluded in Table|2] It should be cautioned that the effective tem- 
peratures of the stars are not well constrained using photomet- 
ric points and should ideally be determined by high resolution 
near-infrared spectroscopy. The derived luminosities thus de- 
pend sensitively on the assumption that the stars are indeed pre- 
main sequence stars. As a consistency check the stellar ages 
and masses have b een derived using the evolutionary tracks of 
ISiess et all (Eo OO). These are also given in Table |2] but are as- 
sociated with large uncertainties except in the cases of SVS 4-9 
and SVS 4-10, which have spectroscopically determined spec- 
tral types. 

In order to convert optical depth and extinction to col- 
umn density, the relation s A^o/th^o = 1.56 x 10 18 cm~ 2 
(Pontoppid an et all2003al) andA^/Ay = 5.6xl0 21 cm^mag -1 
(Vuongetal. 2003) are used. The water ice column den- 
sity conversion uses a band strength of 2 x 10~ 16 cmmolec 1 
dderakines et al.1 19951) . 

3.3. Derivation of the extinction and water ice optical 
depth toward SVS 4-12 

The determination of the extinction and water optical depth to- 
ward SVS 4-12 must be discussed in detail, since this source 
is not detected in the J- and H- band images from UKIRT and 
since the water ice band is heavily saturated. Only a K-bwA 
point, the 3.8-4.2/mi spectral slope and an uncertain 6.7 /urn 
ISOCAM point are available for the extinction determination. 
In principle, the lack of near-infrared photometric points pre- 
vents an accurate determination of extinction if the source 
has a significant infrared excess. However, an extinction of 
Ay = 26 mag gives a good fit to a Rayleigh-Jeans slope of a 
stellar spectrum, without any infrared excess (see Fig.|6j. The 
presence of any excess will result in a smaller extinction, since 
the dereddened SED in that case will be shallower. Therefore, 
the derived extinction is a conservative upper limit. However, 
the simultaneous good fit of the /T-band point and the slope of 
the ISAAC spectrum indicates that SVS 4-12 has a naked pho- 
tospheric SED at least to wavelengths of 4/vm, although the 
current data do not allow us to rule out the possibility that SVS 
4-12 has some minor infrared excess at longer wavelengths. 

The optical depth of the water ice band toward SVS 4-12 is 
determined by scaling the water band from SVS 4-9 to the SVS 
4-12 red wing. In general, all the water ice bands observed to- 
ward the SVS 4 sources have very similar shapes, i.e. one water 
band can generally be scaled to fit well with another at the 10% 
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SVS 4-5 


18:29:57.63 


+01:13:00.2 




0.03 ± 0.01 


2.9 ± 0.3 


130 + 10 


242 + 50 


494 ± 100 


SVS 4-6 


18:29:57.69 


+01:13:04.4 


0.5 ±0.1 


3 ±0.3 


10.3 ± 1.0 


38 + 4 


40 + 6 


84 ±20 


SVS 4-7 


18:29:57.84 


+01:12:27.8 


0.015 ±0.005 


0.6 ±0.1 


4.3 ± 0.4 


15 + 2 


7 + 2 


7±2 


SVS 4-8 


18:29:57.88 


+01:12:37.7 


0.04 ± 0.01 


1.8 ±0.2 


14.4 + 1 


56 + 6 


51 + 10 


56 ± 10 


SVS 4-9 (EC 95) 


18:29:57.92 


+01:12:46.0 


0.21 ±0.05 


9.7 ± 1.0 


68.6 ±7.0 


155 + 15 


147 + 30 


89 + 20 


SVS 4-10 (EC 92) 


18:29:57.88 


+01:12:51.1 


0.6 ±0.1 


9.9 ± 1.0 


41.4 ±4.0 


111 + 11 


152 ± 40 


700 + 70 


SVS 4-12 


18:29:57.22 


+01:12:58.4 






0.2 ± 0.02 


26 + 3 


24 + 5 





All flux densities are given in units of mJy. 



Table 2. Physical parameters for SVS 4 sources 



Source 


Lboi 


Teg 


stellar mass 


Age 




1"HbO 


Nh 2 o 


H2O ice abund; 




[L ] 


[K] 


[M ] 


[Myr] 


[mag] 




[10 18 cm" 2 ] 


xl0~ 5 w.r.t. H 2 


SVS 4-2 


0.7 


4750 


1.5 


17 


4.2 


0.7 ±0.1 


1.1 ±0.2 


9.3 ± 1.4 


SVS 4-3* 


4.9 


4250 


1.0 


0.5 


18.0 


3.5 ±0.5 


5.5 ±0.8 


10.9 ± 1.6 


SVS 4-4 


0.9 


4250 


1.0 


3 


4.9 


0.85 ±0.1 


1.3 ±0.2 


9.7 ± 1.2 


SVS 4-5 


38 


4500 


3.5 


0.3 


21.0 


3.5 ±0.3 


5.5 ± 0.5 


8.6 ±0.8 


SVS 4-6 


0.4 


4750 


0.8 


42 


4.9 


0.7 ±0.1 


1.1 ±0.2 


9.1 ± 1.4 


SVS 4-7 


1.5 


4500 


1.3 


3 


10.5 


2.0 ±0.3 


3.1 ±0.5 


10.6 ± 1.8 


SVS 4-8 


4.6 


4500 


1.4 


0.8 


10.8 


1.7 ±0.2 


2.7 ±0.3 


8.8 ± 1.2 


SVS 4-9 


27 


4750 


3.5 


0.4 


10.5 


1.6 ±0.1 


2.5 ± 0.2 


8.5 ±0.8 


SVS 4-10 


1.2 


3750 


0.5 


0.8 


6.2 


1.1 ±0.1 


1.7 ±0.2 


9.9 ± 1.2 


SVS 4-12* 


28 


5000 


3.5 


0.5 


26 


8± 1 


12.5 ±2.0 


17.2 ±2.4 



Estimated relative uncertainty of 5%. 

* Not detected in J, making the derived spectral type and age uncertain. 



level. The best fit of th 2 o = 8 + 1 for SVS 4-12 using the water 
band of SVS 4-9 is shown in Fig. [7] This value does not change 
significantly if other water bands from SVS 4 are used. Thus 
if the extinction is overestimated due to an infrared excess in 
SVS 4-12, the measured abundance of water ice relative to the 
dust column density will be a strict lower limit. 

4. Relation to the envelope of SMM 4 

The projected distances of the individual SVS 4 stars to the 
central submillimeter position of SMM 4 are 15-40", corre- 
sponding to 3 750 - 10 000 AU at a distance of 250 pc. SMM 4 
is also associated with a large scale molecular outflow oriented 
NW-SE. Consequently, the projected south-eastern lobe of the 
outflow passes through the SVS 4 region. Therefore, there is a 
possibility that the lines of sight toward some SVS 4 stars may 
probe the outer envelope of SMM 4 and/or material which has 
been processed by shocks created by the outflow. 

It is therefore an important issue to determine the depth of 
the SVS 4 cluster relative to the envelope and outflow of SMM 
4. Clearly, much of SVS 4 share the line of sight with the SMM 
4 envelope as seen in Fig. Q] However, it is not clear whether 
the SVS 4 stars are situated in front of, inside, or behind the 



envelope of SMM 4. Only the two last scenarios allow for a 
study of the ices in the envelope of a class protostar since the 
ice absorption bands do not probe what lies behind the individ- 
ual SVS 4 stars. The last scenario would require the presence 
of a second dense cloud in front of SVS 4, which can produce 
the ice absorption and extinction corresponding to H2 column 
densities of up to 7 x 10 22 ctrT 2 for SVS 4-5 and SVS 4-12. The 
Serpens clo ud was mapped in a r ange of optically thin molec- 
ular lines bv lOlmi&Testil J20021 who found a column density 
through the SVS 4 region of 6 x 10 22 ctrr 2 within a large beam 
of ~ 50". However, this column density was measured using 
C 18 0, a molecule known to be partially frozen out at low tem- 
peratures and high densities. The direct measurement of the CO 
ice abundance discussed in Sec. l5.3l shows that the gas-to-solid 
ratio of CO is 2-5 toward SVS 4, assuming a total abundance 
of CO of 2 x lO -4 . However, if the CH 3 OH and C0 2 ices have 
been chemically formed from hydrogenation and oxidation of 
CO molecules, the fraction of CO molecules directly observed 
to be removed from the gas-phase increases by more than 50% 
(see Taking the freeze-out into account, the H2 column 
densities may increase to 12 x 10 22 cirr 3 in the most extreme 
case. The observed column densities can therefore not com- 
pletely rule out the presence of a dense cloud containing SVS 
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Fig. 4. The spectral energy distributions of the SVS 4 stars, which were used to define a continuum for the water bands. The 
extracted water ice bands are shown next to each SED on an optical depth scale. On the left, the dashed lines show the Kurucz 
stellar atmosphere model, the dotted line shows the added infrared excess and the full line is the sum of the two contributions. 



4 and the observed ices, which is unrelated to SMM 4. A high 
resolu tion map (~ 20") of the C ls O (2-1) line bv lWhite et ail 
il 199.4 shows that SVS 4 is located on the edge of the molec- 
ular ridge of the south-eastern clump in Serpens, but shows 
no evidence of a separate, extended clump of dense gas with 
very high column density associated with SVS 4. Additionally, 
the C 18 line has a single gaussian shape with a FWHM of 
~ 2 km s" 1 in the direction of SMM 4, indicating that no other 
large scale velocity component is present. It will be assumed 
that SVS 4 is not associated with a foreground cloud, although 
it should be noted that this issue is not fully resolved. 



To decide whether SVS 4 can be placed behind or signif- 
icantly inside the envelope of SMM 4, a physical model for 
the envelope was calcu lated using the procedure described in 
Jeirg ensen et alJ 1 20021) . The model was fitted to the SCUBA 
850jum map bv bavis et all i 19991) an d the ISO-LWS spec- 
trum extracted bv lLarsson et alJ J2000I) . The resulting model 
has a power law density profile with index -1.6 and an outer 
radius of slightly less than 10000AU where the temperature 
drops to 9 K. The total luminosity is 1 1 L Q , the density at 
1000 AU is 5 x 10 6 cm" 3 and the total mass is 5.8 M . These 
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Fig. 5. Continued from Fig. |4] 



values are largely consi stent with earlier simple models by 
iHogerheiide et alJdl999b . 

The projected extinction through the envelope of SMM 4 
calculated from the model density profile is shown in Fig. [3] 
Clearly, the extinction through the SMM 4 envelope is too high 
for all the SVS 4 stars to lie entirely behind the SMM 4 enve- 
lope. However, SVS 4-5 and SVS 4-12, which both have mea- 
sured extinctions of Aj ~ 25 mag, lie very close to the Aj -26 
contour. This places them effectively behind the envelope if all 
the extinction is caused by SMM 4. Conversely, SVS 4-2 and 
SVS 4-6 must in any case lie in the frontal part of the envelope 
because their measured extinctions are much less than those 
predicted by the SMM 4 model. It is interesting to note that 



this behaviour is expected if SVS 4 is a roughly spherical clus- 
ter directly embedded into the envelope of SMM 4. In other 
words, since the angular extent of the cluster is similar to the 
depth of the envelope (~ 10000AU), it is expected that some 
sources lie in front of the envelope and some behind. This ge- 
ometry is thus consistent with the drastic variation of extinction 
from source to source on very small angular scales. 

5. The distribution of ice toward SVS 4 

5.1. Constructing an ice map 

Once the ice column densities and infrared extinctions are de- 
termined, maps can be constructed showing the spatial distri- 
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Fig. 7. Determination of the optical depth of the water ice band 
toward SVS 4-12. The dots show the parts of the SVS 4-12 
L-band spectrum which are detected with more than 5cr. The 
curves are the spectrum of SVS 4-9 scaled to optical depths of 
7, 8 and 9 to indicate the allowed range of depths of the SVS 
4-12 water ice band. 



bution of these quantities toward SVS 4. Some caveats apply to 
such maps. Since only pencil-beam lines of sight are probed, 
no information about the ice content of the cloud between the 
lines of sight can in principle be inferred. However, under the 
assumption that the observed lines of sight resolve any varia- 
tions in ice abundance, a map can be constructed using interpo- 
lation techniques. In this work, the IDL 2-dimensional kriging 
interpolation algorithm is used. 




40 20 -20 -40 

4-7 aroseo 



Fig. 8. Contours showing the projected extinction in the /-band 
through the model envelope of SMM 4 compared to the posi- 
tions of the SVS 4 lines of sight. The solid lines indicate the 
range of observed position angles of the outflow associated 
with SMM 4. 



The creation of an ice map is similar to the construction 
of extinction map s of dark clouds using background stars (e.g. 
lAlves et alJEoOll) . Due to sensitivity limitations, background 
stars can often not be used to obtain ice maps of high spatial 
resolution. In this article, a densely packed collection of stars 
embedded in the molecular cloud rather than background stars 
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are used, and consequently great care has to be taken in inter- 
preting the resulting distribution of ice. One problem concerns 
the interpretation of differences in column densities, since em- 
bedded stars are likely to probe different depths through the 
cloud. This problem will not affect maps of abundances of the 
ice relative to extinction or other ice species. The second caveat 
concerns the interaction of the embedded stars with the ice in 
their immediate surroundings. We will show in the following 
section that this is unlikely to cause problems for stars of low 
luminosity (L ~ 1 L G ). Conversely, maps can be used to study 
the interaction of individual embedded stars with local ice. 

5.2. Water ice 

5.2.1 . The water ice map 

Maps of the water ice column density, extinction and water ice 
abundance toward SVS 4 are shown in Fig. [5] The water ice 
abundance map is simply the ratio between the water ice col- 
umn density map and the extinction map. It thus shows the line 
of sight averaged abundance of water ice relative to the refrac- 
tory dust component or H2 column density in front of the SVS 
4 stars. It can be seen that the extinction and the column den- 
sity of water ice vary with factors of 5 and 10, respectively, on 
angular scales of only ~ 10". The largest extinction and water 
ice column density are observed toward SVS 4-12. The abun- 
dance of water ice relative to molecular hydrogen is also seen 
to peak sharply toward SVS 4-12 by a factor of 1.9. Within the 
framework of the envelope model of SMM 4, the line of sight 
toward SVS 4-12 passes through a gas density of 4 x 10 5 cirT 3 
and a temperature of 1 1 K. Furthermore, the abundance rises 
slightly in the south-western part of the cluster, parallel to the 
dens e ridge in the south -eastern Serpens core seen in molecular 
gas dWhite et all 1995b . 

The relation between the infrared extinction and the water 
column density is show in Fig. [TO] The derivative of this rela- 
tion basically measures the water ice abundance as a function 
of depth into the cloud. Note the important distinction to the 
absolute line of sight abundance, which is simply the ratio of 
ice column density to H2 column density. The derivative of the 
abundance relation is a measure of the local ice abundance av- 
eraged over the angular dimension of the observed cloud. An 
absolute abundance is an average over local abundances and 
does for example not remove any contribution from bare grains 
along the line of sight. The absolute ice abundance can thus be 
seen as a lower limit to the local ice abundance. 

The lines of sight toward the SVS 4 sources exhibit a linear 
relation or constant water ice abundance between Aj = 5-25 
except for the line of sight toward SVS 4-12, which is over- 
abundant in water ice by a factor of 1.9 relative to the other 
SVS 4 sources. The best-fitting least-squares line for SVS 4 
between Ay =4-25 and excluding SVS 4-12 is: 

t h ,o = (0.07 ± 0.09) + (0.15 + 0.01) x Aj (1) 

The slope of the derived relation corresponds to a local abun- 
dance of water ice relative to H2 of (9 ± 1) x 10~ 5 . Extrapolating 
the linear relation to small extinctions yields a negative inter- 
cept with the A y axis of -0.5 + 0.6. However, the negative in- 




ARCSEC 

Fig. 9. Top: Map of the water ice optical depth toward SVS 4. 
Middle: Map of Ay of the same region. Bottom: Map of the 
abundance of water ice relative to H2. Filled circles indicate 
the positions of the measured sight-lines. SMM 4 is centered 
on (0,0). 
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The corresponding relation for the Taurus m olecular cloud 
ofT H ,o = 0.072 x (A v- 3. 2) is included in Fig.fToll Whittet et al. 



2001). The dense cloud extinction law derived by these authors 
in Taurus, Ay = rEj-x with r = 5.3, was used to convert the 
optical extinction to A y, giving th,o = 0.24 x (Aj + 0.95). The 
conversion is relevant for extinctions above the ice threshold. 
The water ice relation for Taurus was determined using sources 
with Ay < 26, roughly corresponding to Ay < 7.5. It is seen to 
be consistent with the SVS 4 points out to Ay ~ 6, but devi- 
ates significantly from the best-fitting line through the SVS 4 
points. According to the Taurus relation, the water ice abun- 
dance in Taurus is 1.35 x 10~ 4 . In other words, assuming that 
the grain size distribution in SVS 4 and in Taurus are identical, 
the ice mantle volume must be on average 40% smaller in SVS 
4 compared to Taurus in order to create the SVS 4 ice abun- 
dance. However, the abundance of water ice toward SVS 4-12 
is significantly higher than that of the Taurus cloud. 

5.2.2. A jump in the abundance of water ice in the 
SMM 4 envelope? 

The very high abundance of water ice observed toward SVS 
4-12 merits further scrutiny. Fig.^2 snows the observed water 
abundance toward each SVS 4 source as a function of the max- 
imum density (or the minimum distance to SMM 4) probed by 
each line of sight. The densities and distances were calculated 
using the model profile of SMM 4 and the measured extinctions 
of the individual sources to estimate their depth in the SMM 
4 envelope. SVS 4-12 is the source which probes the deepest 
into the envelope at a minimum distance to the center of SMM 
4 of 4700 AU and a maximum density of 4.3 x 10 5 cirT 3 . In 
the figure, a sharp jump in the water ice abundance is evident 
just below 5000 AU, while the abundance remains constant at 
larger distances. However, the presence of the jump depends on 
only one source, SVS 4-12. As discussed in Sec. 13. 31 we regard 
the abundance measurement toward SVS 4-12 relative to those 
of the other SV S 4 sources as robust o r possibly a lower limit. 
As discussed by J0rge nsen et alJ J2004I) . the amount of freeze- 
out is expected to increase significantly at the envelope radius 
where the density is high enough for the freeze-out timescale 
to become shorter than the age of the core. For a typical life- 
time of ~ 10 s yr, this is expected to occur at densities of a 
few xlO 5 cm -3 , comparable to the density at which the jump 
occurs. Of particular interest is the apparent sharpness of the 
jump, which seems to be no more than 500 AU wide. If real, 
this provides strong constraints on such models and justifies 
further ice mapping of protostellar envelopes. Unfortunately, 
the data presented here only give constraints on the abundance 
of water ice toward SVS 4-12 and additional deep mid-infrared 
spectroscopy is required to study the abundance jump in other 
ice species. 

5.2.3. Caveats 

To test the conclusions reached by analysing the th20 - Ay- 
relation, sources of systematic errors should be explored. For 
example, a shallower infrared extinction law would result in a 



Fig. 10. Relation between extinction and water ice column den- 
sity for the SVS 4 sources (filled circles) and the 3 unre- 
lated Serpens YSOs, EC 90A+B and EC 82 (open squares). 
The solid line indicates the best-fitting line through the SVS 
4 points, excluding SVS 4-12. The dashed line is the best- 
fitting line through the three comparison Serpens sources and 
the dotted lin e is the relation measured for the Taurus cloud by 
IWhittet et aJ] ll200ll) (see text). 



tercept is not statistically significant and the relation could also 
pass through the origin. Note that a negative intercept with the 
Ay axis, if significant, translates into the presence of water ice 
without a refractory dust component. This is an unlikely sce- 
nario, and more plausible explanations of a negative threshold 
would be that a systematic error plays a role or that the rela- 
tion is no longer linear for small A y. Significantly, the SVS 4 
sources seem to exclude a positive threshold for the water ice 
in the th,o - A y-relation of Ay ~ 1 such as found in general 



for molecular clouds 



re.g.lWhittetetalJll988h and for Serpens 
: Hodarjc 1989). The two sources SVS 4- 



in particular llEiroa & 
3 and SVS 4-7 have a slightly, but significantly, higher water 
ice abundance compared to the best fitting line of about 20%. 
These two sources are responsible for the higher water abun- 
dance in the south-western part of SVS 4 as seen in Fig. [9] 

For comparison, lines of sight toward three other young 
stellar sources from the south-eastern clump in Serpens (EC 
90A+B and EC82), unrelated to SVS 4, have also been plot- 
ted in Fig. These objects were observed and analysed us- 
ing the same approach as for the SVS 4 sources. It is found 
that the other Serpens sources are systematically and signif- 
icantly under-abundant in water ice compared to SVS 4 and 
result in a relation with a small, but positive intercept with the 
A y axis of Ay = 0.4 +1.3. This is roughly consistent with pre- 
vious determinations of the water ice-extinction relation for 

■ III 
Serpens (Eiroa& Hodarjp 1989), but also consistent with the 

relation passing through the origin. This may be an indication 
of a higher degree of ice evaporation near these sources or a 
lower ice formation efficiency outside the densest parts of the 
Serpens core. 
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Fig. 11. Relation between the water and methanol ice abun- 
dances and the SMM 4 envelope model. The bottom axis indi- 
cates the maximum density encountered by each line of sight. 
The top axis indicates the corresponding minimum distance to 
the center of the SMM 4 envelope. The dashed line sketches 
the sharp jump in water ice abundance around 5000 AU from 
the center of SMM 4. The dotted line sketches the jump in 
methanol ice abundance observed at a radius of ~12000AU 
relative to the center of SMM 4. Outside this radius only strict 
upper limits on the abundance of solid methanol are available 
(see Sec.l5"4l. 

systematic underestimate of the extinction of up to 10%, suffi- 
cient to cause the relation to pass through the origin. However, 
as discussed in Sec. 13.21 a small index for the extinction law 
is not favoured. Also, this effect would influence the three 
comparison sources. The errors in the individual water opti- 
cal depths are often larger than 10%, but systematic effects are 
estimated to be smaller. The dominant systematic sources of 
error in the optical depths are related to the adopted model for 
the intrinsic SED. This error can be estimated by comparing 
the shape of the water ice band for sources with and without 
a large infrared excess. In Figs. @] [5] and [6] it is seen that all 
water ice bands have similar shapes regardless of the shape of 
the background source SED. This is a clear indication that the 
derived optical depth is reasonably model independent. 

Finally, it can be expected that the stars heat the dust in 
their immediate vicinity to temperatures high enough to evapo- 
rate water ice. The question is if this dusty, but ice-free region is 
large enough to significantly affect the observations. A zone of 



dust with no ice mantles around each star would create a thresh- 
old in the th 2 o —A ./-relation. We used a 1 -dimensional radiative 
transfer model of a 50 L Q star embedded in a cl oud of con- 
stant density, calculated using the DUSTY code dlvezic et alJ 
1999), to estimate the threshold created by the radiation of the 
star. These models show that the temperature rises above an as- 
sumed water sublimation temperature of 90 K only within a ra- 
dius of 300 AU. For a star of 1 L , the sublimation temperature 
is reached at 50 AU. The contribution to the extinction from 
bare grains is then Aj 0.7 mag for the most luminous stars in 
S VS 4 and Aj ^ 0.1 for the least luminous stars. Correcting for 
this contribution using the calculated luminosity of each star 
brings the th,o - A y-relation slightly closer to the Taurus rela- 
tion. 

No line of sight is expected to pass through the ice evapo- 
ration zones from any of the other stars. If the density around 
the stars is increasing inwards, i.e., if each star in SVS 4 is as- 
sociated with a remnant envelope, the column density of warm, 
bare grains and thus the ice threshold will increase. Conversely, 
this scenario is unlikely since no significant ice threshold has 
been observed. A full map of CO ice toward the SVS 4 stars 
will be much more sensitive to moderate heating of the enve- 
lope material due to the higher volatility of CO ice. The limited 
CO ice observations toward SVS 4 sources described in Sec. 
15. 31 do indicate that a significant fraction of the dust in the re- 
gion presently has temperatures of ~ 20 K, in agreement with 
moderate heating of the cloud and no significant density gradi- 
ents near the SVS 4 stars. 



5.3. CO ice 

An additional constraint on the temperature of the ice can 
be found from high resolution spectroscopy of CO ice in 
the stretching vibration m ode around 2139 cm -1 = 4.67 //m. 
Laboratory experiments bv lCollings et"aTl ll2003) indicate that a 
layer of pure CO ice will partly migrate into an underlying layer 
of porous water ice upon warmup. The different types of CO ice 
can be distinguished by the s hape s of their corresponding spec- 
tral profiles (Tiel ens et alJll99ll) . IPontoppidan et alJ 112003b) 
suggested that the relative column densities of the pure CO 
ice band at 2139.9 cm 1 and of CO ice trapped in water ice 
at 2136.5 cirT 1 are a sensitive indicator of the thermal history 
as well as current temperature of interstellar ice due to the mi- 
gration effect. In this scenario, a strong band of pure CO ice in- 
dicates the presence of ice with temperatures less than ~ 20 K, 
while a strong band of water-rich CO indicates that a signif- 
icant amount of ice at some point in time has been heated to 
temperatures higher than 20 K. 

Three of the SVS 4 stars have been o bserved in CO ice 
(see Table|3}. The spectrum of SVS 4-10 bv lChiaret all ill 994 
is not of high enough quality to extract quantitative informa- 
tion about column densities of different CO ice components, 
although a rough estimate is given in Table [3] The observed 
total line of sight column densities of CO ice in SVS 4-5 and 
SVS 4-9 are 2.7 x 10 18 cm" 2 and 2.5 x 10 18 cm" 2 , respectively. 
However, since SVS 4-5 has an extinction in the /-band of 22.8 
magnitudes, while SVS 4-9 has an extinction of 10.5, the total 



12 



K. M. Pontoppidan et al.: Mapping ices in protostellar environments on 1000 AU scales 



line of sight CO abundance for SVS 4-9 (8.5 x 1(T 5 ) is twice 
that of SVS 4-5. Interestingly, the line of sight abundances of 
pure CO toward SVS 4-5 and 9 are almost identical at 3 x 10 5 , 
while a factor of almost four distinguishes the abundances of 
CO in water. The estimated colum n density of pure CO toward 
SVS 4-10 from lChiar etail ( 1994) is consistent with a constant 
local abundance of pure CO in SVS 4. 

Although it is not possible to identify any unambiguous re- 
lationship using only two points, some interesting suggestions 
can be made. Clearly, the amount of water-embedded CO ice 
on a grain varies significantly through the cloud. This stands in 
sharp contrast to the pure CO ice, which seems to have a con- 
stant local abundance, based on three sources. This is difficult 
to explain from simple migration of CO from the pure form to 
the water-embedded form upon warmup, but also does not con- 
form to a scenario where the water-embedded CO has formed 
together with the water mantle. In the case of efficient migra- 
tion, a large component of CO in water should be reflected 
by a correspondingly smaller component of pure CO. In the 
simultaneous formation case, the water-embedded CO might 
be expected to show a constant local abundance. Instead, the 
CO ice abundances may point to a more complex thermal his- 
tory of the cloud material. A constant local abundance of pure 
CO might indicate recent freeze-out, while the abundances of 
water-embedded CO could indicate thermal processing of an 
earlier CO layer. Also, a processing scenario other than heat- 
ing which drastically changes the water-rich CO abundance 
could easily change the abundances of other ice species, such 
as methanol and CO2. This is not observed to the same de- 
gree as for CO (see Secs. l5.4l and l5.5> In this picture, the early 
thermal processing might have been caused by the formation 
of the SVS 4 cluster, while the recent freeze-out of CO might 
be related to the formation of SMM 4. Alternatively, the CO 
ice bands may show processing from interaction of the SMM 4 
outflow with the envelope material. A complete map of SVS 4 
in CO ice will be needed to explore this question further. 
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Fig. 12. The four lines of sight in SVS 4 where a spectrum of 
good quality of the 3.53^/m CH3OH band could be extracted. 
All spectra are depicted on the same optical depth scale. The 
spectra of SVS 4-5 and 9 are new medium resolution (R=3300) 
observations; the other two spectra are taken using the low res- 
olution mode of ISAAC yielding a resolving power of R=600. 
The vertical dashed line indicates 3.535 /mi. 

Table 4. CH3OH ice observations toward SVS 4 sources 



Source 


N(CH 3 OH) 


Abundance 




[10 18 cnr 2 ] 


[10~ 5 ] w.r.t. H 2 


SVS 4-5 


1.9 ±0.2 


3.0 ± 0.4 


SVS 4-8 


0.80 ± 0.08 


2.6 ± 0.4 


SVS 4-9 


0.64 ± 0.06 


2.2 ± 0.3 


SVS 4-10 


0.56 ± 0.06 


3.2 ±0.5 



5.4. Methanol ice 
5.4.1 . Abundance 

One of the primary objectives of mapping ices toward SVS 4 
was to explore the large excess of methanol ice known to be 
present in this region. Fig. ^] shows the spectra of the four 
lines of sight from which a high-quality 3.53 fim methanol band 
could be extracted. The remaining lines of sight do not have 
sufficient signal-to-noise ratios to allow a meaningful estimate 
of the column density of solid methanol. In Fig.^]the relation 
between extinction and the column density of methanol ice is 
shown. The best-fitting line to the tch,oh _ A ./-relation is: 

tch 3 oh = 0.06 + 0.1 1 + (0.07 ± 0.01) x Aj. (2) 

The relation gives a local abundance of methanol ice of 
(2.5 + 0.4) x 10~ 5 relative to H2. The measured threshold value 
is consistent with mag and may therefore be identical to the 
water ice threshold. There is a significant scatter in the rela- 
tion, which is not seen for the water ice band for the same four 
sources. The methanol abundances in SVS 4-5 and SVS 4-10, 



both of which are located in the northern part of SVS 4, are 
slightly higher than for SVS 4-8 and SVS 4-9. For comparison, 
strong upper limits on the methanol abundance for EC 90 A 
and B are shown in Fig.[0] In absolute terms, the upper limits 
for these sources are 3 - 4x 10~ 6 . This large difference between 
SVS 4 and EC 90 is particularly interesting since EC 90 is lo- 
cated only 75" north of SVS 4-5. Therefore, the methanol ice 
must reside in a region no more than 2 arcminutes in extent. 
This region contains only SVS 4, and the outflow and envelope 
of SMM 4. There are thus strong indications that the high ef- 
ficiency of methanol formation in the ice mantles in SVS 4 is 
related to the protostellar envelope of SMM 4. 

5.4.2. Methanol band profile 

The detected methanol bands toward SVS 4 sources are indis- 
tinguishable within the noise. In order to maximise the signal- 
to-noise, a weighted average of the four detected methanol 
bands has been calculated. The resulting 3.53/mi methanol 
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Table 3. CO ice observations toward SVS 4 sources 



Source 


N(Pure CO) 


N(CO in water) 


Abundance (pure CO) 


Abundance (CO in water) 


Reference 




[10 I8 cnT 2 ] 


[10 18 cnT 2 ] 


[10~ 5 ] w.r.t. H 2 


[10~ 5 ] w.r.t. H 2 




SVS 4-5 


1.8 ±0.4 


0.9 ±0.1 


2.8 ±0.7 


1.4 ±0.2 


a 


SVS 4-9 


0.9 ±0.15 


1.6 ±0.1 


3.1 ±0.6 


5.4 ± 0.6 


a 


SVS 4-10 


0.4 ± 0.2 


1.0 ±0.5 


3±2 


6±3 


b 



References: a)|Pontoppidan et al. 1 2003a), b)IChiar etall( 1994) 
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Fig. 13. Relation between extinction and methanol ice column 
density for SVS 4. Filled circles indicate SVS 4 sources, while 
open squares are upper limits for the unrelated Serpens YSOs 
EC 90 A+B. The solid line is the best-fitting line through the 
SVS 4 points. 

band profile is shown in Fig. [H] where it is compared to an 
ISAAC spectrum of the well- studied met hanol band toward 
the massive YSO W 33 A (e.g. lDartois etal1ll999t) . While the 
methanol bands from SVS 4 are similar in shape to that of 
W 33A, they almost completely lack the so-called 3.47 yum 
band, which is prominent toward W 33A. A laboratory spec- 
trum of a m ixture of H?0:CQ7 :CH^O H= 1:1:1 ices deposited 
at 10 K from Ehrenfreund et al. ( 1999) is compared to the SVS 
4 methanol band. The main band at 3.53//m is seen to be 
well matched by the laboratory profile in position and width. 
The secondary bands at 3. 32-3. 42 /im are in general harder to 
match with a laboratory spectrum. This may be due to an un- 
certain continuum determination in the blue wing of the band, 
since most of the SVS 4 sources are very faint below 3.35 /im. 
However, even the brightest source, SVS 4-9, shows no sign of 
excess absorption below 3.35/vm. Some substructure absorp- 
tion bands may be seen at 3.40 and 3. 47 /mi, the last of which 
is discussed in more detail in Sec. 15.61 

5.5. C0 2 ice 

Abundant CO2 ice has been observe d toward SVS 4 using 
ISOCAM-CVF dAlexander et al.ll2003l) . Unfortunately, many 
of the SVS 4 sources are blended in the CVF images, espe- 
cially at the longer wavelengths, and accurate decompositions 
are very difficult due to the spatial undersampling of the de- 




O.lO l 1 1 1 I 

3.40 3.50 3.60 3.70 

Wavelength [micron] 



Fig. 14. Weighted average of the observed methanol bands 
from SVS 4-5, 4-8, 4-9 and 4-10 (black dots) compared to an 
ISAAC spectrum of the same region toward the massive YSO 
W 33A (grey dots). The solid line is a laboratory spectrum of 
a H 2 0:C0 2 :CH 3 OH=l:l:l mixture at 10 K. The bands from 
solid H2CO which may be superposed on the main methanol 
band around 3.53 /jm are indicated. 

tector. However, our ISOCAM photometry of the high resolu- 
tion broad-band imaging shows that a few sources dominate 
the flux at 14.3 yum. For instance, SVS 4-5 is 6 times brighter 
than SVS 4-6 with which it is blended. Similarly, SVS 4-10 is 
8 times brighter than SVS 4-9. In principle, the optical depth 
of an absorption band toward two blended sources will depend 
on the ratio of continuum fluxes as well as on the absolute op- 
tical depths of the two lines of sight. In particular, the depths 
of optically thick bands may be significantly underestimated by 
filling-in from another line of sight. The observed optical depth 
of an ice band, r b. s , from two blended sources, S 1 and S 2, is 
given by: 

exp(-r obs ) = expC-nXl + l/R]- 1 + exp(-T 2 )[l + R\~\ (3) 

where R is the ratio of the continuum flux level of S 1 to the 
continuum flux level of S 2. The filling-in effect is most se- 
vere if the two sources are equally bright (R = 1) and one 
source has an optically thick ice band (t > 1). For instance, 
the observed optical depth will never exceed ln(2) 0.7 if the 
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ice band is not present in one of the blended, equally bright 
sources. Therefore, if the CO2 ice bands are not optically thick, 
it should be possible to use the ISOCAM spectra of the two 
bright sources to estimate reasonably accurate optical depths. 

The observed optical depths of the SVS 4-5/6 blend and the 
SVS 4-9/10 blend are 0.5 and 0.25, respectively. In the extreme 
and unlikely scenario that SVS 4-6 and SVS 4-9 have no C0 2 
ice bands, the optical depths of SVS 4-5 and SVS 4-10 would 
be corrected to 0.61 and 0.29, respectively. Conversely, if the 
fainter sources have a very deep ice band, the actual optical 
depths would be corrected to as little as 0.35 and 0.13. Since 
the extinction of SVS 4-5 is 5 times higher than the extinction 
of SVS 4-6, a large correction is expected for this line of sight 
due to filling-in from SVS 4-6. A smaller correction is expected 
for SVS 4-10, since the second source has an extinction which 
is only 1.7 times higher. Assuming that the local CO2 abun- 
dances for the blended sources are identical, the most likely 
C0 2 ice optical depths for SVS 4-5 and SVS 4-10 are then 0.60 
and 0.23, corresponding to absolute abundances of 2.0 x 10" 5 
and 3.5 x 10~ 5 , respectively. Although more uncertain, the CVF 
spectrum of SVS 4-8 gives a CO2 ice abundance of 2 x 10~ 5 . 
There are thus some indications of a slightly varying CO2 ice 
abundance across SVS 4. It is interesting to note that the two 
sources with the lowest absolute abundance of CO2 ice are the 
same sources with the highest absolute abundance of methanol 
ice. Although not conclusive, this may indicate a competition 
between the formation of CO2 and CH3OH from CO through 
oxidation and hydrogenation, respectively. 

5.6. Presence of formaldehyde? 

The abundance of formaldehyde (H2CO) ice in star forming 
regions has been a much debated issue. If the methanol ice 
has been formed through successive hydrogenation of CO, a 
significant abundance of formaldehyde is expected in the ice 
(Watanab e et alJl2003l) . Lines of sight toward massive young 

stars have been found to have fairly small abundances of 

1 . I 

formaldehyde of at most 5% relative to water (Dartois et al. 
Il999l iKeane et alJEoOll) . H 2 CO has several strong bands, but 
they all suffer from blending with other strong features. The 
most isolated band is the V4 mode at 3 .47 /mi. There is a sharp 
band at this position which is common when methanol ice is 
abundant. However, the corresponding v\ band of H2CO at 
3. 54 /mi would create a significant shoulder on the red side 
of the 3.53/mi methanol band, which is not clearly observed 
in SVS 4 (see Fig. I14>. This is the same conclusion reached 

— • 1 j — n 

for other sources by Dartois et al. ( 1999). Formaldehyde ice 
has its strongest bands centered at 5.8 /mi and 6.69 /mi. These 
bands are unfortunately blended with the 6.0 //m water ice band 
and the 6.85 /im band, which has still not been unambiguously 
identified. In Fig.lBlthe ISOCAM CVF spectrum of the 5-8 /mi 
region toward the SVS 4-5/6 blend is shown to explore whether 
a significant abundance of formaldehyde is consistent with the 
shape of the ice bands in this region. The SVS 4-5 line of sight 
has the deepest ice bands as well as the highest signal-to-noise 
in the CVF pointing. The ice bands are compared to laboratory 
spectra of pure water ice and pure formaldehyde ice at 10 K. 



-0.1 
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Fig. 15. ISOCAM CVF spectrum of the 5-8 /mi region of the 
SVS 4-5/6 blend. The curves are laboratory spectra of pure 
water ice (dotted line) and pure H2CO ice (dashed line). Both 
spectra have been convolved to a resolving power of R=35 to 
match the CVF spectrum. The water band has been scaled to 
match the 3.1 /mi band observed with ISAAC toward SVS 4-5. 
The H2CO band has been scaled to the depth of the 3.47 //m 
band. The solid line is the sum of the two laboratory spectra. 

The two laboratory spectra have been convolved to the resolu- 
tion of the CVF spectrum and scaled to optical depths of 3.5 
and 0.1 of the 3.1 //m water band and the 3.47 /mi formalde- 
hyde band, respectively. This corresponds to an abundance of 
formaldehyde ice of 1 x 10" 5 or about 10% relativ e to water ice, 
using a band strength of 0.96 x 10 17 cmmolec -1 (Schutte et al. 
1993) for the 5.8/im band. The CVF resolution significantly 
under-resolves the two formaldehyde bands, so a spectrum of 
higher resolution will significantly enhance the contrast of the 
formaldehyde to the underlying 6.0 and 6.85 //m bands. 

It is found that the water ice band cannot account for the 
entire 6.0 /im band, but falls short by about 50%. This is a well- 
known problem and probably indicates that one or more addi- 
tional carriers contribute to the band. Typically, it is found that 
the 6.0/mi band is twice as deep as it should be when com- 
paring to the 3.1 /mi band dOibb & Whittetll2002t) . The reason 
that a slightly smaller discrepancy is found here may be due 
to veiling by SVS 4-6. According to Eq. [5] a blended source 
with a shallow ice band a factor of 10 fainter than the primary 
source will change the optical depth from an observed 0.37 to 
an intrinsic 0.42. The presence of a formaldehyde component 
indeed gives a better match to the blue wing of the 6.0 /zm band. 
However, the 6.69/im band does not seem consistent with the 
shape of the 6.85 /zm band. We conclude that while the spectra 
toward the SVS 4 stars are not consistent with a high abundance 
of formaldehyde, an abundance of ~ 5% relative to water ice is 
possible. Higher resolution spectroscopy in the 5-8 /mi region 
is required to convincingly detect formaldehyde ice. 

5.7. Interaction of the outflow with the ice 

SMM 4 is associated with a large bipolar outflow. Estimates 
of the position angle of the outflow range between 180° 



K. M. Pontoppidan et al.: Mapping ices in protostellar environments on 1000 AU scales 



15 



Table 5. Summary of ice abundances in the outer envelope of 
SMM4 



Species 


Abundance 


Abundance 




w.r.t. H 2 


xl0~ 5 w.r.t. H 2 


H 2 


1 


9-17 


CH 3 OH 


0.28 


3 


Pure CO 


0.3 


3 


CO in water 


0.1-0.6 


1-5 


C0 2 


0.2 


2.0 - 3.5 


H 2 CO 


£0.05 


< 0.5 


XCN 


< 0.003 


< 0.02" 



" IPontoppidan et alj f2003b) 

llHogerheiide et al.lll999h and 135° ddarav et aljEoOol The 
south-eastern outflow lobe is red-shifted, while the north- 
western lobe is blue-shifted. An outflow is expected to evap- 
orate at least part of the ice mantles through shock-induced 
sputtering. 

Observations of rotational lines of gaseous methanol have 
shown that the gas-phase abundance of methanol in the SE lobe 
is enhanced by a factor of 40 relative to typical abundances 
in quiescent molecular gas of CO/CH30H= 3 x 1CT 5 , while 
the a bundanc e in t he NW lobe is enhanced by a factor of 330 
JGarav et alJ EoOO). It is common to see such enhanced abun- 
dances in the gas phase of methanol toward bipolar outflow 
sources. The most favourite interpretation is that the methanol 
originates in icy mantles, which have been shock-evaporated 
by the outflow. However, it is important to note that the peak 
column density of gas-phase methanol toward the SE lobe of 
SMM 4 is only 1.4 x 10 15 cirT 2 , which translates to a gas-to- 
solid ratio for methanol of ~ 2 x 10~ 3 . Consequently, even in 
lines of sight with strong gas-phase enhancement, the vast ma- 
jority of the methanol molecules remain as ice in the grain man- 
tles. It is shown in Sec. 15.41 that the abundance of methanol ice 
is enhanced by an order of magnitude in the SMM 4 outflow re- 
gion. If the grain mantles in the specific case of the SMM 4 out- 
flow were to have a typical abundance of methanol of less than 
a few percent relative to water ice rather than the observed 25%, 
the line of sight gas-phase enhancement of methanol would be 
less than a factor of 4. In this case, it is doubtful that the outflow 
would have been recognised as having an enhanced gas-phase 
methanol abundance. Extrapolating this result to other proto- 
stellar outflows with strongly enhanced gas-phase methanol 
abundances could indicate that the ice-phase methanol abun- 
dance is commonly enhanced in protostellar envelopes. In most 
chemistry models, the ice mantles are assumed to form during 
the pre- and protostellar stages, either by direct freeze-out of 
molecules from the gas (e.g., CO) or grain-surface che mistry 
of acc reted species (e.g., H2O from accreted O). tBergin et alJ 
1998) proposed an alternative scheme in which ice mantles 
form in the wake of low-velocity shocks due to the outflow im- 
pacting on the envelope. At the high shock temperature, most of 
the gas-phase oxygen is driven into H2O, resulting in very high 
H2O ice abundances if all of this H2O freezes out before it is 
chemically altered into other species in the cold post-shock gas. 



Such a scenario could be tested by comparing high-resolution 
maps of gas-phase H2O with ice maps of stars behind outflow 
lobes such as presented here. 

5.8. Absolute ice abundances and gas-phase 
depletion 

The abundances relative to H2 of the different ice species ob- 
served in the outer envelope of SMM 4 are summarised in Table 
[5] The species constitute all of the most abundant products of 
surface oxygen chemistry according to current ice mantle mod- 
els. It is therefore of interest to consider the abundances in re- 
lation to the gas phase components of the protostellar envelope 
and surrounding dense cloud. For instance, since the abundance 
of chemical products believed to originate in CO frozen out 
from the gas is very high, a considerable amount of the CO will 
be bound in the form of other molecules once the ice is returned 
into the gas phase upon warm-up. This will to some extent in- 
validate gas phase abundances determined relative to CO even 
in warm gas. In other words, contrary to the gas phase, CO is a 
chemically active molecule in the solid phase. Toward SVS 4, 
the abundance of molecules formed from CO (CO2, CH3OH, 
H2CO) is up to 7x 10~ 5 . This is up to 2 times as much as the CO 
ice remaining in the grain mantles. Therefore, in some cases, 
only 1/3 of the CO originally frozen out is returned to the gas 
phase as CO. The total abundance of frozen-out CO molecules, 
including those which have hydrogenated or oxidised to other 
molecules is up to 1.4 x 10~ 4 relative to H2, corresponding to a 
high depletion from the gas-phase. 

Another interesting number to mention is the total abun- 
dance of oxygen observed in the ice. For SVS 4, it is as much as 
2.7 x 10~ 4 with respect to H2 and perhaps as much as 3.5 x 10 4 
toward SVS 4-12. In comparison, the gas phase abundance of 
atomic oxygen in the diffuse interstellar medium within 5 00 p c 
of the Sun is O/H = (3.19 + 0.14) x 10" 4 dMever et alJll998h . 
Thus, about half of the available oxygen is bound in ice in the 
dense parts of the Serpens core. A similar conclusion holds for 
carbon. The ice abundance and composition may therefore be 
expected to vary considerably in regions with significantly dif- 
ferent abundances of heavy elements. One example may be in 
star-forming regions in the Large Magellanic Cloud, which is 
known to have a sub-solar abundance of heavy elements. 

6. Conclusions 

We have presented VLT-ISAAC L-band spectroscopy of most 
of the stars in SVS 4, a dense cluster of low-mass YSOs lo- 
cated near the class protostar, SMM 4. These observations 
have been coupled with archival data from ISO and UKIRT to 
produce a detailed view on the distribution of ices in the imme- 
diate environment of a very young protostar. The main conclu- 
sions are: 

- It is shown that SVS 4 is likely located inside the south- 
eastern part of the outer envelope of SMM 4. The clus- 
ter members are distributed evenly throughout the envelope 
such that the lines of sight toward some stars probe the en- 
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tire depth of the envelope. This is in particular true for SVS 
4-12. 

- A 30" x 45" map of the water ice abundance with a spa- 
tial resolution of 6" has been constructed. It shows that 
water ice is distributed throughout the SVS 4 cluster with 
an abundance of 9 + 1 x 10~ 5 relative to molecular hydro- 
gen. The abundance of water ice rises sharply to at least 
1 .7 x 1 (T 4 toward SVS 4- 1 2, which is the line of sight prob- 
ing closest to the center of SMM 4. The ice along this line 
of sight accounts for up to 50% of the available oxygen 
in the interstellar medium. This is a clear indication of an 
enhanced efficiency of ice formation in the inner parts of 
protostellar envelopes. 

- Along with the water ice, methanol ice is present in SVS 4 
with a fairly constant high local abundance of 2.5 + 0.4 x 
10~ 5 or 28% relative to water ice. The methanol ice is con- 
fined to SVS 4, and the upper limit to the methanol ice 
abundance just 75" = 19 000AU away from the center of 
SVS 4 is < 3 x 10~ 6 or < 5% relative to water ice. This 
testifies to a strongly enhanced efficiency of the formation 
of solid methanol in the envelope of SMM 4. 

- Other ice species show normal abundances relative to wa- 
ter ice. There is, however, tentative evidence for a varying 
abundance of CO2 ice across SVS 4. 

- Mapping of ices at a spatial resolution comparable to that 
of gas-phase emission observations is now possible with 
current instrumentation on 8 m class telescopes. Due to the 
fortunate alignment of a dense cluster of young stars with 
a protostellar envelope it was possible to achieve the high 
spatial resolution of 6". However, a spatial resolution of 
30 - 60" is possible in the more general case. 

- Maps of the abundances of H2O, CH3OH and CO ices rel- 
ative to the column density of refractory dust or H2 are po- 
tentially powerful probes of the physical history of molec- 
ular cores. Further ice mapping is essential to explore the 
possibilities of using ice as a physical probe. Additionally, 
it is important to couple maps of ices with corresponding 
gas-phase maps and models of the density and temperature 
of the dust to obtain a complete picture of the chemistry 
in particular of the saturated species. Indeed, about half of 
all molecules apart from H2 not bound in refractory dust is 
in the form of a chemically active ice for most of the life- 
time of a dense core. The Spitzer Space Telescope allows 
mapping at higher sensitivities in ice species not observable 
from the ground, such as CO2 and CH4. 
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